The shell model approach 
to the rotating turbulence 



M. Reshetnyak"' 6 '*, B. Steffen c 

Institute of the Physics of the Earth, Russian Acad. Sci, 123995 Moscow, Russia 0, 

Research Computing Center of Moscow State University, 
119899, Moscow, Russia 

Central Institute for Applied Mathematics (ZAM) 
of Forshungszentrum Jiiulich, Germany 



Abstract 

Applications of the shell model of turbulence to the case of rapidly rotating bodies 
are considered. Starting from the classical GOY model we introduce the Coriolis 
force and obtain a ~ k~ 2 spectrum for 3D hydrodynamical turbulence for the free 
decay regime as well for the regime with external forcing. Additional modifications 
of the GOY model providing a realistic form of the helicity are proposed. 
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1 Introduction 



The influence of rotation on the properties of the hydrodynamical turbulence 
is of the great importance. This problem appears in the various geophysical 
and astrophysical applications [l( and it requires special treatment. So far the 
direct 3D numerical simulations, where implementation of the Coriolis force 
is trivial, cannotprovide a long enough inertial range of spectrum to reveal 



its scaling laws |2], thus different approaches are needed. 

Though the possibility of transition from 3D-turbulence state to the two- 
dimensional turbulence due to rotation was already predicted by Batchelor 
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many years ago l3j], only now is the qualitative description of this processes 
developing |^|, [5|. Having in mind the similar situation in the MHD tur- 
bulence p], where the magnetic field plays the same role as the rotation in 
reducing the dimension of the problem, the authors proposed that the Coriolis 
force introduces a new characteristic time. This time should be used instead of 
the characteristic turn-over time based on the Kolmogorov's estimate. Then, 
instead of Kolmogorov's —5/3 slope the spectrum will have a —2 slope. It ap- 
pears that this approach is in agreement with the direct numerical calculations 
and the experiments 0. The formal simplicity of this phenomenological 
approach attracts us to implement it in the more complicated models of tur- 
bulence and to check the results for self-consistence. For this aim we use the 
well-known homogeneous isotropical GOY shell model (see overview in [9]) 
and modify it to the case of rotation. The proposed model is tested for the 
regimes of the free-decay turbulence and for the regime with external forcing. 
We also consider a situation where a non-zero average helicity is generated. 
This regime finds its application in the mean-field dynamo problems. 



2 Basic equations 



The evolution of conductive incompressible liquid (V • V = 0), rotating with 
the angular velocity Q, can be described by the Navier-Stokes equation: 

/<9V \ 
i? l— + (V-V)Vj =-VP + F + £V 2 V. (1) 



Eq.(l) is scaled with the large scale of the body L; velocity V, pressure P 
and time t are measured in f]/L, prf / L 2 and L 2 /ri, where 7] is a magnetic 
diffusion, p is density, R a = r]/2QL 2 and E = u/2QL 2 are the Rossby and 
Ekman numbers, respectively and v is the kinematic viscosity. (This kind of 
scaling takes its origin in geodyamo problems 0.) The force F = F c + f 
includes the Coriolis effect: 

F c = -l z x V (2) 



and a prescribed part f . Here l z is the unit vector along the axis of rotation. 

It is worth to note that, in the limit of the vanishing viscosity and absence of 
the external force (f = 0), equation (1) conserves kinetic energy E^ = V 2 /2. 
For the two-dimensional flow the other integral of motion is the enstrophy 
Eq = (V x V) 2 and for 3D case the hydrodynamical helicity l? = V-(Vx V) 
is conserved (see for details, e.g., 
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3 Shell model equations 



The idea of the shell model approach is to mimic the original equation (1) 
by a dynamical system with n max complex variables U\,U2, ■ ■ ■ ,u nmax which 
represent the typical magnitude of the velocity on a certain length scale. The 
Fourier space is divided into n max shells and each shell k n consists of a set of 
wave vectors k such that ko2 n < |k| < fco2 n+1 , and thus each shell corresponds 
to an octave of the wave numbers. Variable u n is the velocity difference over 
length ~ l/k n , so that there is one degree of freedom per shell. The coupling 
between shells is considered to preserve the main symmetries and properties of 
equation (1). Of course, the shell model is a simplified form of the description of 
the turbulence, nevertheless, it is a reasonable tool in the studies of turbulence. 
Here and after we will refer to the GOY model |9||, which has the form: 



R° d dt ~ Roikn (y , n+l U n+2 2 M n-l U n+l ^4 ' u n-2 u n-l) ,„n 

-Ek 2 n u n + F n , 



where e is a free parameter, * is the complex conjugate. 

In the inviscid limit (E — ► 0) and free forcing (F = 0) equation (3) has two 
integrals. The first integral is the kinetic energy 

= £ Kl 2 . (4) 



The other integral is the so-called generalized enstrophy 

F = £[sgn(e-l)]X {e) Kr, (5) 

n 

with a(e) = — log 2 (|e — l|/2). For e > 1 H is always positive 

H = H = Y,nk^\u n \ 2 . (6) 



For e = 5/4 its dimension coincides with dimension of the enstrophy, which is 
an integral of motion in 2D problem. That is the reason why e = 5/4 usually is 
associated with the two-dimensional space. Cases with e < 1 correspond closer 
to the helicity integral and can have different sign. In the case of e = 1/2 the 
dimension of H is equal to the dimension of the hydrodynamical helicity, 
so this case corresponds to 3D space. It is known that the 3D shell model 
can reproduce Kolmogorov's spectrum ~ k~ 5 ^ 3 . The situation in 2D space 
(e = 5/4) is more complicated. Application of an external force at some wave 
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number kf gives rise to two spectral regimes: for k > kf the direct spectrum for 
the enstrophy, and the spectrum for the inverse energy cascade with the slope 
similar to Kolmogorov's —5/3 for k < kf. The slope of the inverse cascade 
spectrum for enstrophy depends on e and can change its value from —3 for 
e = 5/4 to —5/3 for e « 10 llj. Here and after we will consider only the direct 
cascades in the 3D case with e = 1/2. 

To introduce effect of rotation in the GOY model we propose that the Coriolis 
force can be written in the form: 

F c = -iC r u, (7) 



where the constant real coefficient C r is introduced for convenience. It is easy 
to see that the work of the Coriolis force (7) is zero {u* n F cn + u n F* n = 0) and no 
additional energy is introduced into the system (3) at any scale. This property 
of the force (7) coincides with the property of the Coriolis force in 3D space 

(!)• 

Having in mind that in derivation of the shell model equations (3) all external 
potential forces, as well as pressure, were already excluded using condition of 
incompressibility (V ■ V = 0), F c corresponds to the curl part of the Coriolis 
force. We will return to this point when the case with the several forces will 
be considered in section 5. 

In the following sections we will consider the applications of the shell model (3, 
7) to the free decay turbulence (section 4) and to the regime with the external 
forcing (section 5). Special attention will be paid to the helicity generation 
produced in such models (section 6). 



4 Free decay hydrodynamic turbulence 



To analize behaviour of the hydrodynamical turbulence without forcing we 
start from the free decay regime (/ = 0). After some intermediate regime for 
the case without rotation (F c = 0) the Kolmogorov's spectrum (—5/3) recov- 
ers, see figures 1, 2. The sharp breaks in the spectra at the large wave numbers 
correspond to the Kolmogorov's diffusion scale k&. The estimate based on the 
balance of the inertial R k n u\ and diffusion Ek^u n terms leads to k^ ~ ^fu n . 

As was predicted in 0, 0], introduction of rotation (F c ^ 0) gives rise to a new 
time scale t& ~ Q" 1 . For R <C 1, r<j is already shorter then the characteristic 
time in the Kolmogorov's regime r^ n ~ k~ 2 / 3 . A simple dimensional analysis 
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leads to the estimate of the rate of the energy dissipation 
e ~ r(k)k A E 2 (k), 



(8) 



where r is the characteristic time and E(k) is the spectral energy density In 
the case of the Kolmogorov's turbulence r ~ (k 3 Ey 1/2 and 

E(k) ~ e 2/3 k~ 5/3 . (9) 



If the effect of rotation is sufficient, then substitution of r ~ Q 1 into (8) leads 
to the rotation spectrum law jlj]: 

E(k) ~ (fie) 1/2 Ar 2 . (10) 



Starting simulations from the initial field obtained in the non-rotating regime, 
after a short time period the original Kolmogorov's spectrum splits into two 
parts with two different slopes, see Fig. 1, 2. The change in the slope (% G 

[10, 191 and k n = 2 n ) corresponds to k n = where C n = 1.22 1.87 0. 

This estimate can be obtained from balancing the inertial term and the Coriolis 
force. If for the large scales the Coriolis term is larger then the non-linear term, 
the spectrum decays as ~ k~ 2 . In this case non-linear term do not depend on 
k and the Coriolis term decays as ~ k^ 1 ^ 2 . The further behaviour depends 
on how long the spectrum is and where the Kolmogorov's wave number k^ 
lies. If kn > k d , then the whole spectrum decays as ~ k~ 2 . In the other 
case (kn < k d ) the Kolmogorov's spectrum —5/3 for k > k n reappears. This 
situation is demonstrated in the figures 1, 2. As the whole kinetic energy of 
the system decays in time, the k^ moves in to the region of the large wave 
numbers. Note that the situation can be complicated by the high order effects, 
as the decay rates of the both spectra are different. 



5 Convection with external forcing 



To consider the evolution of the system for time periods longer than the char- 
acteristic decay time one needs to provide some source of energy to the system 
or to introduce an external force /. Similar to the case with the Coriolis force, 
/ in (3) corresponds to the curl part of the force. 

The results of simulations of the system (3) over the time period t = 10 3 with 
a prescribed force / = 10~ 2 (1 + i) at ko = 1 without rotation are presented in 
Fig. 3. Starting from an arbitrary initial velocity field, the system (3) comes 
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to the statistically stable state with Kolmogorov's energy distribution, similar 
to the free decay case for the moment where the energy is comparable. 

As was already mentioned above, the direct introduction of rotation setting 
C r = 1 in (7) is contradictive to the physics of the problem and to the original 
equation (1). It is easy to see that for the regime of the fast rotation (R a <C 
1) in the original Navier-Stokes equation (1), when V < R , the balance 
between the pressure and the potential parts of the forces / and F c holds. 
Such a balance of the pressure and the Coriolis force, which takes place, e.g., 
in the Earth's liquid core (E = 10~ 15 , R Q = 4 • 10~ 7 ), is called the geostrophic 
balance (in the present dimensionless units regime the Earth's core situation 
corresponds to R Q V < 10~ 3 ). Exclusion of the pressure requires exclusion of 
the potential parts of all the forces, too. Then, the remaining curl parts of 
the forces are already of the same order as the inertial term and it is these 
parts which are in the r. h. s. of (3). These considerations can be formulated 
as follows: 

{C rn -R fc n | tl n | , -Rofenjltnl <C \Uji\ (11) 



In other words, condition (11) means that for all wave numbers where the 
Coriolis force is larger then the non-linear term it must be compensated by 
the pressure, and its curl part (F c ) must be of the same order as the non-linear 
term. (An example of violation of the condition (11) will be also considered 
below.) The suggestion (11) is also supported by the recent paper 12||, where 
the scaling law ~ k~ 2 for the geostrophycal balance was derived. 

The results of simulations with rotation and C r defined by (11) are presented 
in Fig. 3. As in the free decay case we observe two regimes with slopes ~ k~ 2 
for the small wave numbers and the Kolmogorov's regime ~ k~ 5 ^ 3 for the 
wave numbers k > k^. It is easy to see that in this model the non-linear term 
has a white spectrum for k < k^. Due to the condition (11), the curl part 
of Coriolis term F c has the same distribution in the wave space. This state 
of equipartition was observed after averaging. Analysis of the phases of the 
non-linear term and the Coriolis term reveals the existence of anticorrelation. 
This is the reason why the spectrum of the rotating turbulence decreases faster 
then Kolmogorov's one. In fact, the Coriolis force partially locks the energy 
transfer from large to small scales. Moreover, the Coriolis force blocks the 
applied external force, so that F c0 = — f c o an d suppresses injection of the 
energy into the system (3, 7, 11). 

To demonstrate this phenomenon we present calculations with C r ten times 
larger then was predicted by (11) (dotted line in Fig. 3). This regime shows 
very strong blocking. Our calculations show that the further increase of the 
Coriolis term leads to a degeneration of the spectrum at all. The solution then 
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has a singular spectrum and the whole energy is concentrated at the wave 
number fc where the force / is applied. 

We also present results of calculations with a prescribed force of the form 
/o = 0.1(1 + i)/uq which allows us to define the exact amount of energy 
injected into the system. It appears that the introduction of rotation leads 
to an increase of the total energy of the system due to a decrease of the 
energy dissipation in the cascade, see Fig. 4. Note that the similar effect in 
3D numerical simulations 0, where the external force was applied at some 
intermediate scale, was caused by the inverse cascade. 



6 Helicity generation 

As follows from (5) in 3D case (e = 1/2) the explicit form for helicity is 

H = Y J {-^TK\u n \ 2 . (12) 

n 

However, the direct application of formula (12) for the estimation of helicity or 
the calculation of a-effect can lead to misleading results. Additional assump- 
tions on the physics of the process are required. To understand the nature of 
the problem we will recall some basics of the mean-field dynamo theory, where 
the models of helicity for the rotating bodies were extensively studied. 

According to 0], the mean helicity in the rotating body is defined as H =< 
V ■ (V x V) >, provided that the gradient of the density exists (< ... > 
means mean average). This problem reflects the statistical inequality of the 
left-side and right-side rotating curls. In its turn, such inequality is caused 
by existence of the selected direction (z), which is related to the global rota- 
tion of the considered body. In absence of rotation the generated helicity has 
stochastic temporal behavior and zero mean level. Moreover, in the case of the 
Kolmogorov's turbulence with E{k) ~ k~ 5 ^ 3 , the helicity H has an increasing 
spectrum ~ fc 1//3 and its amplitude will be defined by the velocity field near the 
Kolmogorov's cut-off wave number 1 . The temporal behaviour of the he- 
licity for the case of the non-rotating homogeneous turbulence (with all other 
parameters are equal to the case with forcing in section 5 and summation in 
(12) is over the whole range of wave numbers) is presented in Fig. 5a. We 
observe high frequency oscillations with the negligible mean level (H ~ 10 2 ). 



1 Note that for the dynamo applications the maximal wave number in the sum (12) 
is limited by the condition r m > 1, where r m is the micro-scale magnetic Reynolds 
number. 
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This case can hardly describe the helicity observed in astrophysical and geo- 
physical bodies with a stable sign over time periods much longer then the 
characteristic convective times |14]. 



Application of the Coriolis force in the form (7) does not change the situation 
essentially. In Fig. 5b we present calculations of H for the previously con- 
sidered regime with rotation (Fig. 1). Even now, when the helicity spectrum 
is white (~ k°) and correspondingly contributions of the high harmonics are 
smaller then in the Kolmogorov's turbulence, the helicity behavior exhibits 
strong oscillations with asymptotically zero mean level (H ~ 10 1 ). 

To estimate which part of helicity is produced in the rotational part of the 
spectrum {k < k^) we limit the cut-off wave number in (12) to k M Ax = k^, 
see Fig. 5c. This filtering decreases the dispersion of the helicity, but still does 
not change its mean level (H ~ 10 1 ). 

To overcome the problem of zero mean level of helicity we have to introduce 
an additional asymmetry in the system (3,7). So far in our homogeneous in- 
compresible shell model the gradient of the density is absent, so we have to 
include the missing effect in the Coriolis force: 

F c = -zC r u(l + d(-l) n ), (13) 

where C\ is a real constant. The results of simulations with C\ = 4, kuAx = k n 
are presented in Fig. 5d. In this case H ~ 140 and is comparable with the 
amplitude of the oscillations. The corresponding spectrum for this regime has 
the same behavior as in the case when the Coriolis force was defined by (7), 
see Fig. 3. 



7 Conclusions 

The proposed shell model approach description to the rotating hydrodynami- 
cal turbulence demonstrates the principal possibility to reproduce the scaling 
law ~ k~ 2 predicted in pj, [^|. This deviation from the homogeneous Kol- 
mogorov's turbulence represents the reduction of the 3D space problem toward 
2D space, where the direct cascade of enstrophy already has ~ k~ 3 slope. In 
terms of the shell models one can solve the Cauchy problem for evaluating 
rotating turbulence and mimic its dynamical behavior in time. As was men- 
tioned, the problem of rotating turbulence has applications to the liquid core 
of th Earth. The estimate of the Reynolds number based on the west drift 
velocity and molecular viscosity gives Re ~ 10 9 , which is out of reach of any 
modern supercomputers in 3D direct numerical simulations. The situation in 
many other astrophysical bodies is even more dramatic where the estimates 
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of the Reynolds number give Re ^> 10 10 [14]. This Reynolds number is no 
problem in the shell model, where the wide range of scales can be covered by 
a few harmonics. That is why the proposed application of the shell model to 
the fast rotation regime is a very promising step in its development. We also 
realize that the mechanism of helicity generation proposed by is not the 
only one and some forms of helicity generation different from (13) are pos- 
sible. The present approach was dictated by its popularity in the mean-field 
dynamo model of the a-effect. We plan to use the shell model developed for the 
computation of the subgrid processes in the 3D large-scale thermal convection 
problem similar to (l5| . 
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Fig. 1. Spectra of the free decay hydro dynamical turbulence, E = 10 , R = 10 . 
Solid line corresponds to the regime without rotation, C T = 0, and line with circles 
to the regime with rotation, C r = 1. 




Fig. 2. Spectra of the free decay rotating turbulence, E = 10 , R Q = 4 • 10 




Fig. 3. Spectra of the turbulence with external forcing /o = 10~ 2 (l + i), E = 1CT 10 , 
R = 1CT 3 . Solid line corresponds to the regime without rotation, C r = 0; line 
with circles corresponds to the regime with rotation, C T = 1 and the Coriolis force 
defined by (7); the dotted line to C r = 10. The line with squares corresponds to the 
modified Coriolis force (13). 




Fig. 4. Evolution of the kinetic energy for the non-rotating (1) and rotating (2) 
turbulence, f = 10 _1 (1 + i)/u* , E = 1(T 10 , R Q = 1(T 3 . 



8000-1 



ffi 



-8000 - 1 
2000-1 

X 




a 



a o 



ffi 0- 



Fig. 5. Time evolution of helicity H. The curve (a) corresponds to the non-rotating 
(Kolmogorov's) turbulence; curve (b) is the regime with the Coriolis force defined 
by (7); case (c) is the same to (b) but with cut-off at kn (only k < kn in (12)) and 
case (d) with the Coriolis force defined by (13) and k < kn. 



